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Abstract. Existing stellar X-ray surveys suggest major problems in our understand- 
ing of the evolution of stellar magnetic activity in solar and late-type stars, reaching 
conflicting conclusions over the rate of decay of X-ray activity and the spectral types 
responsible. We are confronting these discrepancies with a new model of the expected 
stellar X-ray luminosity distribution, combining a Galactic population synthesis model 
with current theories for rotational spin- down and the rotation - activity relation for the 
stellar magnetic dynamo. Here we test our model using new observations of the stellar 
content of the Chandra COSMOS survey, for which 60 new stellar X-ray sources are 
identified from the thin disk and Galactic halo populations. Our model is in approxi- 
mate agreement with the observed X-ray luminosity distribution and the distribution of 
spectral types responsible. However, slight differences in the form of the X-ray lumi- 
nosity distribution exist that may hint at problems in our understanding of stellar X-ray 
emission. 



1. Introduction 

Solar and late-type stars emit X-rays through a magnetically confined plasma, or corona, 
at several million Kelvin. The corona is believed to be powered by t he stellar dynamo . 



which itself is powered by differential rotation within the star (e.g. Skuma nichI ll 9721 : 



IPallavicini et ak 19 81; No yes et al. 1984) . The observed decrease in stellar X-ray lu- 
minosity of several orders of magnitude be t ween the zero age main sequence (e.g. 
Feiaelson et al."2002'; 'Flacc omio etaP l2003l : IWright et al.ll2010bh and solar age (e.g. 



Peres et al. 20 00) has therefore been attributed to the rotational spin-down of the star. 



though a consistent picture has yet to emerge. X-ray observations of stellar clusters 
up to a billion years in age have been u sed to stu dy the age - rotation and rotation - 
activity relationships (e.g. Barnes.2003. ; ,Pizzolato et al.,,2003.) . However, the dispersal 
of stellar clusters as they age means that very few stellar clusters older than a billion 
years exist, and many of these are faint and distant. Age estimates for individual main- 
sequence field stars are difficult to derive independently and therefore the evolution of 
both rotation and stellar activity for older stars is poorly understood. 

Attempts to overcome this probler n by using large sa mples of field stais have of- 
ten resulted in conflicting conclusions. iGiidel et akl d 19971) studied a sample of nearby 
solar-type stars age d 1 - 1 Gy r and found that the X-ray luminosities decayed as 



Lx t while Micela (200 21) could f i nd no evidence for a clear decay law over 
a similar age range, and Feigel soii et akl (l200 4l) estimated a decay law of Lx oc t~^ 



from a sample of faint high Galactic latitude main-sequence stars. Other studies have 
attempted to understand the evolut ion of stellar X-ra y emission by comparing predic- 
tions from models (e.g. XCOUNT. lFavata et ak|[T992l) with X-ray surveys of field stars. 
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However, these methods have also uncovered discrepancies such as an observed ex - 
cess of yellow (G-K dwarf) stars compared to model predictions ( Favata et al.lll988h . 
and evidence for an exc ess of young main- sequence stars suggesting either a recently 
high star forrnation rate ( Micela et al.ll 19931) or uncertainties in the binary star fraction 
(iMicela et alJl2007h . 

These discrepancies are no doubt partly influenced by the small but diverse range 
of stellar X -ray samples used. The se include wide-field surveys from Einstein and 
ROSAT (e.g. ISchmitt & Lie&3l2004b that may be biased toward bright and nearby thin- 
disk stars as well as deep surveys with the Ch andra and XMM-New ton observatories 
that include many distant or faint sources (e.g. iFeigelson et al.ll2004t) . Between these 
two extremes d eep and wid e-field surveys such as the Chandra Cosmic Evolution Sur- 
vey (COSMOS. IeMs et al.ll2009) off'er the balance necessary to produce large samples 
of stellar X-ray sources that are also deep enough to uncover intrinsically faint sources. 

In this contribution we discuss a new model for predicting the X-ray emitting stel- 
lar population in any field of view of the Galaxy. We test this model against a new and 
well-constrained catalog of stellar X-ray sources detected in the Chandra COSMOS 
survey. 



2. XStar: a stellar X-ray emission model 

To further our understanding of the evolution of stellar dynamo activity we wish to 
compare the X-ray luminosity distribution of a well constrained sample with that pre- 
dicted from current models of the X-ray evolution of late-type stars. Doing this for 
individual stars is difficult because of uncertainties in d eriving accurate stellar ages. An 
alternative method that has been used in the past (e.g. Favata e t"anil992l: iMicela et al 



[l993) is to compare the observed X-ray luminosity distribution with predictions from 
a model that takes into account the stellar populations along the line of sight and the 
spectral type dependent age - rotation - activity evolution of late-type stars. 

To do this we have a compiled a new X-ray population synthesis model, XStar, 
that combines the latest theories of Galactic structure, rotational evolution and stellar 
dynamo activity to make predictions about the X-ray emitting stellar population in any 
field of view. The ingredients of this model are as follows: 

• A model of the spatial distribution and properties of stars in the Galaxy along 
a partic ular line of sight. For this we use the TRILEGAL population synthesis 
model (iGirardi et al.ll2005h . which can be calibrated using deep optical observa- 



tions of the field in question. 

A theory for the spectral-type dependent stellar spin-down as stars age. For this 
we use the fr amework outlined bv Barnes (2003) and in particular the rotational 
isochrones of Meibom et al. I (l2009h . 



The X-ray activity - rotation relationship quantified by iPizzolato et all (2003), 



using a new empirical determination of the convective turnover time based on 
rotation period measurements of field and cluster stars. 

The product of this model is a probability distribution function of the X-ray emit- 
ting stars along the line-of-sight, including spectral types, distances, ages, rotation peri- 
ods, and X-ray luminosities. This distribution can then be compared to the observed dis- 
tribution of these quantities to allow the components making up the model to be tested. 
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Where discrepancies exist between model and observations, different quantities can 
be compared to test different aspects of the model such as the Galactic star formation 
history (through the distribution of distances) or the binary star fraction (through the 
spectral type distribution). A fully calibrated model can then be used to make predic- 
tions of the X-ray emitting foreground population that may contaminate studies of star 
forming regions (e.g. .Wright & Drake 20 09) or young clusters ( Pillitteri et ah 2006) , or 
be used to assess the late-type stellar cont ribution to the extended X-ray emission from 
nearby galaxies (e.g. iBoroson et al.ll201(il) . 



3. Testing the model on the Chandra COSMOS field 

To test the XStai" model we use the new catalog of ste llar X-ray sources in the Chandra 



COSMOS survey presented by Wri ght et all (I2010al) . The deep and abundant optical 



and near-IR photometry for the field allows the Galactic starcou nts model to b e accu- 
rately calibrated, while the completeness fractions estimated by iPuccetti et al.l (l2009h 
can be used to modulate the resulting X-ray distribution functions. 

3.1. New stellar X-ray sources in the Chandra COSMOS field 

The Cosmic Evolution Survey (COSMOS. IScoville et al.ll2007bl ) is a deep and wide ex- 
tragalactic survey designed to probe the medium redshift galaxy and active galactic nu- 
clei populations. The field h as been observed at nearly all wavelengths (e.g. lCapak et al. 



20071 : IScoville et al.ll2007ah pr oviding a large m ulti-wavelength catalog. The Chandra 



survey of the COSMOS field (lElvis et al.ll2009h has imaged an area of -0.9 deg2 with 



a uniform exposure of -160 ks over the centr al 0.5 deg^ and ~80 ks over the outer re- 
gion. A detailed source detection procedure ( Puccetti et al.ll2009 h resulted in a catalog 



of 1761 sources with well defined completeness fractions as a function of both X-ray 
band and survey area. 

Optical and near-IR identifications for 1750 of these sources were made by'Civano| 

(l2010h . including a number of stellar sources. AA^right et al. (2010a) used photometry 



and spectroscopy to confirm the stellar nature of 60 of these sources, obtaining spectral 
classifications, photometric distances and extracting X-ray properties such as plasma 
temperatures and luminosities. The sources uncovered span the full spectral sequence 
of coronal X-ray emitters, from late F-type to mid M-type, with distances ranging from 
30 pc to 12 kpc, extending the known coverage of the Lx - distance plane to greater 
distances and higher X-ray luminosities than ever before (see Figure [T|l. 



The majority of sources are more luminous than the contemporary Sun (iPeres et al. 



2000l) and the spectral-type dependence of the X-ray luminosities is in good agreement 



with previous studies (e.g. ICovey et all 120081) . The most distant of the sources are 



likely to be the most luminous members of the X-ray - emitting Galactic halo popula- 
tion, which appears to comprise both low-activity spectrally hard sources, that may be 
emitting through thermal bremsstrahlung, as well as a number of highly active sources 
in close binaries. Two sources are identified at distances of ~10 kpc, which also have 
the highest X-ray luminosities in our sample, Lx ~ 10^^~^^ erg s"^ We consider the 
possibility that these are either mis-identifications, confused sources or background 
fluctuations but conclude that they are more likely to be bona fide stellar X-ray sources. 
Their high X-ray luminosities, despite their considerable age, suggests that they are in 
close binaries that have been kept active through tidal interactions and the tapping of 
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Figure 1. Lx distance for this and ot her samples dHiinsch et al] 119991; 
Schmitt & LiefS 120041: iFeigelson etaP 12004 lAgiieros et alJ '2009': 'Cove v et alJ 
2008h . Dashed lines show estimated sensitivity limits IFeigelson et al. 2005|) for 
the three Chandra surveys (in the same color as their symbols). Also show are the 
range of distances (above) and Lx (right) spanned by each sample (in the same color 
as their symbols), with the mean and standard deviations marked as dots and dashes 
respectively. 




Figure 2. Starcounts for the COSMOS field (black) and the best-fit TRILEGAL 
model (green) in the SDSS and 2MASS bands, with completeness limits and Icr 
Poisson uncertainties marked. 



orbital angular momentum to sustain strong dynamo activity. See lWright et all COlOah 



for a full analysis and discussion of this sample and its X-ray properties. 
3.2. Calibration of the Galactic starcounts model 

A catalog of all stellar sour ces in the COSMOS field was compiled us ing data from the 
COSMOS optical catalogs ("Caoak et al."2 007l : lMcCracken et allEm d). the Sloan Digi- 



tal Sky Survey (SPSS , York et al. 2000) and the Two Micron All Sky Survey (2MASS 



Skrutskie et al.l 20061) . All the photometry was converted onto a single photometric 



system based on the SDSS and 2MASS systems, and then filtered for galaxies using 
morphological and spectral energy distribution cuts based o n the known l ocus of the 



main-sequence in the various two-color diagrams analyzed by ICovey et al.l ( 20071) . The 



product of these cuts was to reduce the original COSMOS catalog of 1 .7 million sources 
initially to ~500,000 point-like sources, and then to only ~75,000 sour ces with stellar 



colors. The TRILEGAL Galactic star counts model (iGirardi et al.ll2005T) was then used 



to fit a simple Galactic structure model to this sample (see Figure |2ll, altering a small 
number of parameters to reflect local differences in Galactic structure parameters that 
produce a better fit to the observations. The result of this population synthesis model is 
a list of stellar sources in the field of view with important properties such as their age, 
spectral type and binarity that can then used as the input in our stellar X-ray emission 
code. 

3.3. Results 

The results of our model are shown in Figure |3]compared to the observed cumulative X- 
ray flux distribution of sources in the COSMOS field. The model provides a reasonably 
good fit to the brightest sources, but under-predicts the number of sources with fluxes 
between 2x10"'^ and 2x 10~'^ erg s~' cm~^. Fainter than this the model begins to over- 
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Figure 3. Cumulative X-ray flux distribution for the COSMOS stellar sample and 
the results of our model (with the COSMOS completeness limits applied). 



predict the number o f sources. How ever, below fx = 10^'^ erg s ' cm"^ the COSMOS 
survey is incomplete dPuccetti et al.ll2 009). and therefore differences between the model 
and observations could equally be attributed to uncertainties in the completeness curves 
for the COSMOS observations (which were determined based on the spectral energy 
distributions of galaxies and may therefore be inaccurate if used for stars). 

Considering only the cumulative flux distribution above the completeness limit, 
the model under-predicts the number of sources at the 2cr level. The spectral type 
distribution and distance distribution of the model and observations are in better agree- 
ment, suggesting that the discrepancy does not lie with some unresolved component of 
the model (e.g. a recently enhanced star formation rate), but with the age - rotation - 
activity relationships. Given that the age - rotation relationship is better calibrated than 
the rotation - activity relationship, the source of the discrepancy is more likely to be the 
latter. A steeper rotation - activity relationship could partly resolve this discrepancy, 
such that in the simple approximation Lx oc \P,^^^, p > 2 (where ;6 - 2 is the standard 
exponent used in the rotation - activity relationship, e.g. lPizzolato et aL 2003h. 

A similar steepening of this relationship was found by iFeigelson et all (|2004|) 
from their smaller sample of Chandra Deep Field North stars. They found a decay 
of X-ray activity with time of Lx ^ t~^, which for slowly-rotating stars spinning down 
Skumanich-style (Vrot «: r^^^) reduces to a rotation - activity relation of Lx «: vj^^f. We 
find that this significantly over-predicts t he number of obser ved stellar X-ray sources. 
An intermediate level of decay found by iGtidel et all (1 19971) was Lx oc r^-^, which is 
equivalent to Lx k v^^, for slow rotators, and ofi'ers a better fit, but still over-predicts 
the number of sources. An intermediate value of /3 ~ 2.5 produces the best fit, though 
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at the moment we can only rule out the standard case of p = 2 with a confidence of Icr. 
Larger samples are clearly necessary to resolve this. 

Whether this steepening of the r otation - activity relati on is universal or is an 
age-related effect is hard to determine. iFeigelson et al.l (120041) consider the possibility 
that magnetic field generation becomes more efficient as stars age or that the scaling of 
coronal densities behaves differently for older stars. However, the fact that this apparent 
steepening of the rotation - activity relation has been observed in a number of samples 
that cover a range of ages (ou r stars are more X-ray luminous and therefore likely to be 
younger than thos e found by IFeigelson et al.ll2004h could indicate that this steepening 
is universal. Mon tesinos et al. "(12001) argue that current mode ls of two-layer interface 
dynamos (e.g. Par ked 19931 : ICharbonneau & M acGregoil ll997l ) do not scale simply as 
v^gf, but are a more complex function of a number of parameters, not all of which 
can currently be measured or determined from models. Empirical estimates of this 
scaling can therefore better inform stellar dynamo models, ultimately helping us better 
understand our own Sun. 



4. Conclusions 

To improve out understanding of the decay of stellar dynamo activity we have compiled 
a new model of stellar X-ray activity, XStar, which incorporates a population synthesis 
Galactic model, the theories of spectral-type dependent rotational spin-down and the 
rotation - activity relationship. The results of the model can be compared to samples 
of stellar X-ray sources from a range of observatories. Here we test the model us- 
ing new observations of 60 stellar X-ray sources identified in the Chandra COSMOS 
field. A comparison of the distribution of X-ray fluxes, spectral types, and distances for 
the observed and modeled populations shows a generally good agreement, but slightly 
under-predicts the number of observed sources. Differences between the two X-ray Hux 
distributions suggests that the our current formulation of the rotation - activity relation- 
ship may not be correct, and a steeper relationship may better describe the observations. 
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